* Wealth index calculation. Lady Ortiz, 08/01/2023. Jordan DHS 2023.
cd 'C:\JO81\Wealth '.
include file='wealth_assets.sps'.
execute.
set tvars=both.

* Macro definitions to perform some recoding automatically.
* Note: Do not use an asterisk in front of a macro name to try to comment it out - it will run any way because of the way the macros are expanded.
* Delete the line instead or change the spelling of the macro in the comment.
* Do not use the names of macros in comments - they will also be expanded too, and likely will fail.

*{Create binary variables based on condition and give label }.
define dichotomize( !positional !enclose('(',')') / !positional !tokens(1) / !positional !tokens(1) )
compute !2=0.
if (!1) !2=1.
variable labels !2 !3.
value labels !2 0 'No' 1 'Yes'.
formats !2 (f1.0).
!enddefine.

*(Separate shared toilets from non-shared toilets).
define shared_fac( !positional !tokens(1) / !positional !tokens(1) / !positional !tokens(1) ).
compute !2 = 0.
do if (QH110=1).
+ if (!1=1) !2 = 1.
+ compute !1 = 0.
end if.
variable labels !2 !3.
value labels !2 0 'No' 1 'Yes'.
formats !2 (f1.0).
!enddefine.

*{Reset missing values to 'does not have', change 2 code to 0}.
define no2zero( !positional !tokens(1) )
if (missing(!1) | !1<>1) !1=0.
value labels !1 0 'No' 1 'Yes'.
!enddefine.

*{Construct Variables}.

* Type of housing unit.
dichotomize (QH100A=1) QH100A_1 "Type of housing unit: Apartment".
dichotomize (QH100A=2) QH100A_2 "Type of housing unit: Dar".
dichotomize (QH100A=3) QH100A_3 "Type of housing unit: Villa".
dichotomize (QH100A=4) QH100A_4 "Type of housing unit: Hut/barrack".
dichotomize (QH100A=6) QH100A_6 "Type of housing unit: Other".


* Source of drinking water.
dichotomize (QH101=11) QH101_11 "Source of drinking water: Piped into housing unit".
dichotomize (QH101=12) QH101_12 "Source of drinking water: Piped to yard/plot".
dichotomize (QH101=21) QH101_21 "Source of drinking water: Spring".
dichotomize (QH101=31) QH101_31 "Source of drinking water: Rainwater".
dichotomize (QH101=41) QH101_41 "Source of drinking water: Tanker truck".
dichotomize (QH101=51) QH101_51 "Source of drinking water: Bottled water".
*ichotomize (QH101=96) QH101_96 "Source of drinking water: Other".


* Type of toilet facility.
dichotomize (QH109=11) QH109_11 "Type of toilet facility: Flush to piped sewer system".
dichotomize (QH109=12 or QH109=13) QH109_12 "Type of toilet facility: Flush to pit latrine/ Flush to somewhere else".
*dichotomize (QH109=13) QH109_13 "Type of toilet facility: Flush to somewhere else".
dichotomize (QH109=21) QH109_21 "Type of toilet facility: Ventilated improved pit latrine".
dichotomize (QH109=22) QH109_22 "Type of toilet facility: Pit latrine with slab".
dichotomize (QH109=23) QH109_23 "Type of toilet facility: Pit latrine without slab/open pit".
dichotomize (QH109=61) QH109_61 "Type of toilet facility: No facility/bush/field".
*dichotomize (QH109=96) QH109_96 "Type of toilet facility: Other".

shared_fac QH109_11 QH109_11_sh "Type of toilet facility: Flush to piped sewer system - shared".
shared_fac QH109_12 QH109_12_sh "Type of toilet facility: Flush to pit latrine - shared".
* hared_fac QH109_13 QH109_13_sh "Type of toilet facility: Flush to somewhere else - shared".
shared_fac QH109_21 QH109_21_sh "Type of toilet facility: Ventilated improved pit latrine - shared".
shared_fac QH109_22 QH109_22_sh "Type of toilet facility: Pit latrine with slab - shared".
shared_fac QH109_23 QH109_23_sh "Type of toilet facility: Pit latrine without slab/open pit - shared".
*hared_fac QH109_96 QH109_96_sh "Type of toilet facility: Other - shared".


* Type of cookstove.
dichotomize (QH117=1) QH117_1 "Type of cookstove: Electricity".
dichotomize (QH117=2) QH117_2 "Type of cookstove: Natural gas".
dichotomize (QH117=3) QH117_3 "Type of cookstove: Kerosene".
dichotomize (QH117=4) QH117_4 "Type of cookstove: Coal, wood".
dichotomize (QH117=95) QH117_95 "Type of cookstove: No food cooked in household".
*ichotomize (QH117=96) QH117_96 "Type of cookstove: Other".

* Household has an independent bathroom.
no2zero QH122A.

* Heat source for home.
dichotomize (QH123=1) QH123_1 "Heat source for home: Central heating".
dichotomize (QH123=2) QH123_2 "Heat source for home: Kerosine or solar heater".
dichotomize (QH123=3) QH123_3 "Heat source for home: Gas heater".
dichotomize (QH123=4) QH123_4 "Heat source for home: Air conditioner".
dichotomize (QH123=5) QH123_5 "Heat source for home: Electric heater".
dichotomize (QH123=6) QH123_6 "Heat source for home: Firewood/charcoal/peat".
dichotomize (QH123=7) QH123_7 "Heat source for home: No need for heating".
dichotomize (QH123=95) QH123_95 "Heat source for home: No space heating in household".
*ichotomize (QH123=96) QH123_96 "Heat source for home: Other".


* Type of fuel for home heat.
dichotomize (QH125=1) QH125_1 "Type of fuel for home heat: Electricity".
dichotomize (QH125=3) QH125_3 "Type of fuel for home heat: Solar air".
dichotomize (QH125=4) QH125_4 "Type of fuel for home heat: Cooking gas".
*ichotomize (QH125=6) QH125_6 "Type of fuel for home heat: Alcohol/ethanol".
dichotomize (QH125=7 or QH125=6 ) QH125_7 "Type of fuel for home heat: Diesel/ Alcohol/ethanol".
dichotomize (QH125=8) QH125_8 "Type of fuel for home heat: Kerosene/paraffin".
dichotomize (QH125=10) QH125_10 "Type of fuel for home heat: Charcoal".
dichotomize (QH125=11) QH125_11 "Type of fuel for home heat: Wood".
dichotomize (QH125=12 or QH125=13) QH125_12 "Type of fuel for home heat: Straw/shrubs/grass/ Agricultural crop".
*ichotomize (QH125=13) QH125_13 "Type of fuel for home heat: Agricultural crop".
dichotomize (QH125=14) QH125_14 "Type of fuel for home heat: Animal dung/waste".
dichotomize (QH125=16 or QH125=96) QH125_16 "Type of fuel for home heat: Garbage/plastic/Other".
dichotomize (QH125=18) QH125_18 "Type of fuel for home heat: Peat".
*ichotomize (QH125=96) QH125_96 "Type of fuel for home heat: Other".


* Type of light at home.
dichotomize (QH126=1) QH126_1 "Type of light at home: Electricity".
dichotomize (QH126=2) QH126_2 "Type of light at home: Solar lantern".
dichotomize (QH126=3) QH126_3 "Type of light at home: Rechargeable flashlight, torch, or lantern".
dichotomize (QH126=4) QH126_4 "Type of light at home: Battery powered flashlight, torch or lantern".
dichotomize (QH126=5) QH126_5 "Type of light at home: Biogas lamp".
*ichotomize (QH126=6) QH126_6 "Type of light at home: Gasoline lamp".
dichotomize (QH126=7 or QH126=6) QH126_7 "Type of light at home: Kerosene or paraffin lamp/ Gasoline lamp".
*ichotomize (QH126=8) QH126_8 "Type of light at home: Charcoal".
dichotomize (QH126=9 or QH126=8) QH126_9 "Type of light at home: Wood/ Charcoal".
*ichotomize (QH126=10) QH126_10 "Type of light at home: Straw/shrubs/grass".
*ichotomize (QH126=11) QH126_11 "Type of light at home: Agricultural crop".
*ichotomize (QH126=12) QH126_12 "Type of light at home: Animal dung/waste".
*ichotomize (QH126=13) QH126_13 "Type of light at home: Oil lamp".
*ichotomize (QH126=14) QH126_14 "Type of light at home: Straw/shrubs/grass".
dichotomize (QH126=95) QH126_95 "Type of light at home: No lighting in household".
dichotomize (QH126=96 or QH126=10 or QH126=14) QH126_96 "Type of light at home: Straw/shrubs/grass/ Straw/shrubs/grass /Other".

* Radio or tape recorder.
no2zero QH132A.
* Television.
no2zero QH132B.
* Satellite.
no2zero QH132C.
* Land telephone.
no2zero QH132D.
* Refrigerator.
no2zero QH132E.
* Freezer.
no2zero QH132F.
* Washing machine.
no2zero QH132G.
* Dish washer.
no2zero QH132H.
* Solar heater.
no2zero QH132I.
* Air conditioner.
no2zero QH132J.
* Fan.
no2zero QH132K.
* Water cooler.
no2zero QH132L.
* Microwave.
no2zero QH132M.
* Digital camera.
no2zero QH132N.
* Watch.
no2zero QH133A.
* Bicycle.
no2zero QH133C.
* Motorcycle or scooter.
no2zero QH133D.
* Animal-drawn cart.
no2zero QH133E.
* Tablet.
no2zero QH133F.
* Computer.
no2zero QH133G.
* Household has internet.
no2zero QH133AA.

* Main floor material.
dichotomize (QH152=11) QH152_11 "Main floor material: Earth/sand".
dichotomize (QH152=31) QH152_31 "Main floor material: Parquet or polished wood".
dichotomize (QH152=32) QH152_32 "Main floor material: Tiles".
dichotomize (QH152=33) QH152_33 "Main floor material: Marble/ceramic".
dichotomize (QH152=34) QH152_34 "Main floor material: Cement".
dichotomize (QH152=96) QH152_96 "Main floor material: Other".


* Main roof material.
dichotomize (QH153=21) QH153_21 "Main roof material: Mud bricks".
dichotomize (QH153=22) QH153_22 "Main roof material: Mud bricks with stones".
dichotomize (QH153=31) QH153_31 "Main roof material: Concrete".
dichotomize (QH153=32) QH153_32 "Main roof material: Zinc".
dichotomize (QH153=96) QH153_96 "Main roof material: Other".


* Main wall material.
dichotomize (QH154=21) QH154_21 "Main wall material: Mud bricks".
dichotomize (QH154=22) QH154_22 "Main wall material: Mud bricks with stones".
dichotomize (QH154=31) QH154_31 "Main wall material: Cement bricks".
dichotomize (QH154=32) QH154_32 "Main wall material: Cut stones".
dichotomize (QH154=33) QH154_33 "Main wall material: Cut stones and concrete".
dichotomize (QH154=34) QH154_34 "Main wall material: Concrete".
dichotomize (QH154=35) QH154_35 "Main wall material: Zinc".
dichotomize (QH154=96) QH154_96 "Main wall material: Other".

*{Members per sleeping room}.
if (hhusual=0) hhusual=hhslept.
if (QH127>0) memsleep=trunc(hhusual/QH127).
if (QH127=0) memsleep=hhusual.
if (missing(QH127) or QH127>=99 or memsleep>=98) memsleep=99.
variable labels memsleep 'Number of members per sleeping room'.
value labels memsleep 0 'Less than 1 per room'.
formats memsleep (f2.0).
missing values memsleep (99).

* Compute urban and rural variables coded (1/0) for filters later.
compute urban=(QHTYPE = 1).
compute rural=(QHTYPE = 2).
variable labels urban 'Urban' / rural 'Rural'.
value labels urban 1 'Urban' / rural 1 'Rural'.
formats urban rural (f1.0).

execute.

* Check on indicator variable creation.
frequencies variables=QHTYPE to land.

* Toilet facility by shared/not shared.
crosstabs /tables=QH109 by QH110.

* Land area by units - if there are separate units - need to convert them to one unit.

* Ladys Ortiz.
* Variables with 0 values .
* QH101_96, qh109_96, qh109_96_sh, QH117_96, QH123_96, QH126_11, QH126_12, QH126_13.

* Variables grouped duw small number of cases 0-4.
*  QH109_13, QH109_13_sh, QH125_6, QH125_13, QH125_96 .
* QH126_6, QH126_8, QH126_10, QH126_14.
  .
frequencies variables=
  QH100A_1 QH100A_2 QH100A_3 QH100A_4 QH100A_6
  QH101_11 QH101_12 QH101_21 QH101_31 QH101_41 QH101_51 
  QH109_11 QH109_12 QH109_21 QH109_22 QH109_23 QH109_61 
  QH109_11_sh QH109_12_sh QH109_21_sh QH109_22_sh QH109_23_sh 
  QH117_1 QH117_2 QH117_3 QH117_4 QH117_95 
  QH122A 
  QH123_1 QH123_2 QH123_3 QH123_4 QH123_5 QH123_6 QH123_7 QH123_95 
  QH125_1 QH125_3 QH125_4 QH125_7 QH125_8 QH125_10 QH125_11 QH125_12 QH125_14 QH125_16 QH125_18
  QH126_1 QH126_2 QH126_3 QH126_4 QH126_5 QH126_7 QH126_9 QH126_95 QH126_96
  QH131A
  QH132A QH132B QH132C QH132D QH132E QH132F QH132G QH132H QH132I QH132J QH132K QH132L QH132M QH132N 
  QH133A QH133C QH133D QH133E QH133F QH133G
  QH132AA
  QH133AA 
  MOBPHONE CHECKACC CREDITCARD EWALLET
  QH152_11 QH152_31 QH152_32 QH152_33 QH152_34 QH152_96
  QH153_21 QH153_22 QH153_31 QH153_32 QH153_96
  QH154_21 QH154_22 QH154_31 QH154_32 QH154_33 QH154_34 QH154_35 QH154_96
  HOUSE LAND memsleep
.

* Turn off weights before all factor analysis.
weight off.

* Name the dataset window for the HH data for use later.
dataset name assets.

save outfile='assets.sav'.


*** Test Factor Analysis.
filter off.
execute.

factor /variables
  QH100A_1 QH100A_2 QH100A_3 QH100A_4 QH100A_6
  QH101_11 QH101_12 QH101_21 QH101_31 QH101_41 QH101_51 
  QH109_11 QH109_12 QH109_21 QH109_22 QH109_23 QH109_61 
  QH109_11_sh QH109_12_sh QH109_21_sh QH109_22_sh QH109_23_sh 
  QH117_1 QH117_2 QH117_3 QH117_4 QH117_95 
  QH122A 
  QH123_1 QH123_2 QH123_3 QH123_4 QH123_5 QH123_6 QH123_7 QH123_95 
  QH125_1 QH125_3 QH125_4 QH125_7 QH125_8 QH125_10 QH125_11 QH125_12 QH125_14 QH125_16 QH125_18
  QH126_1 QH126_2 QH126_3 QH126_4 QH126_5 QH126_7 QH126_9 QH126_95 QH126_96
  QH131A 
  QH132A QH132B QH132C QH132D QH132E QH132F QH132G QH132H QH132I QH132J QH132K QH132L QH132M QH132N 
  QH133A QH133C QH133D QH133E QH133F QH133G
  QH132AA
  QH133AA 
  MOBPHONE CHECKACC CREDITCARD EWALLET
  QH152_11 QH152_31 QH152_32 QH152_33 QH152_34 QH152_96
  QH153_21 QH153_22 QH153_31 QH153_32 QH153_96
  QH154_21 QH154_22 QH154_31 QH154_32 QH154_33 QH154_34 QH154_35 QH154_96
  HOUSE LAND memsleep
  /missing meansub
  /analysis
  QH100A_1 QH100A_2 QH100A_3 QH100A_4 QH100A_6
  QH101_11 QH101_12 QH101_21 QH101_31 QH101_41 QH101_51 
  QH109_11 QH109_12 QH109_21 QH109_22 QH109_23 QH109_61 
  QH109_11_sh QH109_12_sh QH109_21_sh QH109_22_sh QH109_23_sh 
  QH117_1 QH117_2 QH117_3 QH117_4 QH117_95 
  QH122A 
  QH123_1 QH123_2 QH123_3 QH123_4 QH123_5 QH123_6 QH123_7 QH123_95 
  QH125_1 QH125_3 QH125_4 QH125_7 QH125_8 QH125_10 QH125_11 QH125_12 QH125_14 QH125_16 QH125_18
  QH126_1 QH126_2 QH126_3 QH126_4 QH126_5 QH126_7 QH126_9 QH126_95 QH126_96
  QH131A 
  QH132A QH132B QH132C QH132D QH132E QH132F QH132G QH132H QH132I QH132J QH132K QH132L QH132M QH132N 
  QH133A QH133C QH133D QH133E QH133F QH133G
  QH132AA
  QH133AA 
  MOBPHONE CHECKACC CREDITCARD EWALLET
  QH152_11 QH152_31 QH152_32 QH152_33 QH152_34 QH152_96
  QH153_21 QH153_22 QH153_31 QH153_32 QH153_96
  QH154_21 QH154_22 QH154_31 QH154_32 QH154_33 QH154_34 QH154_35 QH154_96
  HOUSE LAND memsleep
  /print univariate initial extraction fscore
  /criteria factors(1) iterate(25)
  /extraction pc
  /rotation norotate
  /method=correlation.


*** Common Factor analysis ***.
** Removed area-specific variables **.
** Agricultural animal variables excluded
** Land area excluded

*** Common Factor Analysis.
filter off.
execute.

factor /variables
  QH100A_1 QH100A_2 QH100A_3 QH100A_4 QH100A_6
  QH101_11 QH101_12 QH101_21 QH101_31 QH101_41 QH101_51 
  QH109_11 QH109_12 QH109_21 QH109_22 QH109_23 QH109_61 
  QH109_11_sh QH109_12_sh QH109_21_sh QH109_22_sh QH109_23_sh 
  QH117_1 QH117_2 QH117_3 QH117_4 QH117_95 
  QH122A 
  QH123_1 QH123_2 QH123_3 QH123_4 QH123_5 QH123_6 QH123_7 QH123_95 
  QH125_1 QH125_3 QH125_4 QH125_7 QH125_8 QH125_10 QH125_11 QH125_12 QH125_14 QH125_16 QH125_18
  QH126_1 QH126_2 QH126_3 QH126_4 QH126_5 QH126_7 QH126_9 QH126_95 QH126_96
  QH131A
  QH132A QH132B QH132C QH132D QH132E QH132F QH132G QH132H QH132I QH132J QH132K QH132L QH132M QH132N 
  QH133A QH133C QH133D QH133E QH133F QH133G
  QH132AA
  QH133AA 
  MOBPHONE CHECKACC CREDITCARD EWALLET
  QH152_11 QH152_31 QH152_32 QH152_33 QH152_34 QH152_96
  QH153_21 QH153_22 QH153_31 QH153_32 QH153_96
  QH154_21 QH154_22 QH154_31 QH154_32 QH154_33 QH154_34 QH154_35 QH154_96
  HOUSE LAND memsleep
  /missing meansub
  /analysis
  QH100A_1 QH100A_2 QH100A_3 QH100A_4 QH100A_6
  QH101_11 QH101_12 QH101_21 QH101_31 QH101_41 QH101_51 
  QH109_11 QH109_12 QH109_21 QH109_22 QH109_23 QH109_61 
  QH109_11_sh QH109_12_sh QH109_21_sh QH109_22_sh QH109_23_sh 
  QH117_1 QH117_2 QH117_3 QH117_4 QH117_95 
  QH122A 
  QH123_1 QH123_2 QH123_3 QH123_4 QH123_5 QH123_6 QH123_7 QH123_95 
  QH125_1 QH125_3 QH125_4 QH125_7 QH125_8 QH125_10 QH125_11 QH125_12 QH125_14 QH125_16 QH125_18
  QH126_1 QH126_2 QH126_3 QH126_4 QH126_5 QH126_7 QH126_9 QH126_95 QH126_96
  QH131A
  QH132A QH132B QH132C QH132D QH132E QH132F QH132G QH132H QH132I QH132J QH132K QH132L QH132M QH132N 
  QH133A QH133C QH133D QH133E QH133F QH133G
  QH132AA
  QH133AA 
  MOBPHONE CHECKACC CREDITCARD EWALLET
  QH152_11 QH152_31 QH152_32 QH152_33 QH152_34 QH152_96
  QH153_21 QH153_22 QH153_31 QH153_32 QH153_96
  QH154_21 QH154_22 QH154_31 QH154_32 QH154_33 QH154_34 QH154_35 QH154_96
  HOUSE LAND memsleep
  /print univariate initial extraction fscore
  /criteria factors(1) iterate(25)
  /extraction pc
  /rotation norotate
  /save reg(ALL com)
  /method=correlation.


*** Urban Factor Analysis.
filter off.
filter by urban.
execute.

factor /variables
  QH100A_1 QH100A_2 QH100A_3 QH100A_4 QH100A_6
  QH101_11 QH101_12 QH101_21 QH101_31 QH101_41 QH101_51 
  QH109_11 QH109_12 QH109_21 QH109_22 QH109_23 QH109_61 
  QH109_11_sh QH109_12_sh QH109_21_sh QH109_22_sh QH109_23_sh 
  QH117_1 QH117_2 QH117_3 QH117_4 QH117_95 
  QH122A 
  QH123_1 QH123_2 QH123_3 QH123_4 QH123_5 QH123_6 QH123_7 QH123_95 
  QH125_1 QH125_3 QH125_4 QH125_7 QH125_8 QH125_10 QH125_11 QH125_12 QH125_14 QH125_16 QH125_18
  QH126_1 QH126_2 QH126_3 QH126_4 QH126_5 QH126_7 QH126_9 QH126_95 QH126_96
  QH131A
  QH132A QH132B QH132C QH132D QH132E QH132F QH132G QH132H QH132I QH132J QH132K QH132L QH132M QH132N 
  QH133A QH133C QH133D QH133E QH133F QH133G
  QH132AA
  QH133AA 
  MOBPHONE CHECKACC CREDITCARD EWALLET
  QH152_11 QH152_31 QH152_32 QH152_33 QH152_34 QH152_96
  QH153_21 QH153_22 QH153_31 QH153_32 QH153_96
  QH154_21 QH154_22 QH154_31 QH154_32 QH154_33 QH154_34 QH154_35 QH154_96
  HOUSE LAND memsleep
  /missing meansub
  /analysis
  QH100A_1 QH100A_2 QH100A_3 QH100A_4 QH100A_6
  QH101_11 QH101_12 QH101_21 QH101_31 QH101_41 QH101_51 
  QH109_11 QH109_12 QH109_21 QH109_22 QH109_23 QH109_61 
  QH109_11_sh QH109_12_sh QH109_21_sh QH109_22_sh QH109_23_sh 
  QH117_1 QH117_2 QH117_3 QH117_4 QH117_95 
  QH122A 
  QH123_1 QH123_2 QH123_3 QH123_4 QH123_5 QH123_6 QH123_7 QH123_95 
  QH125_1 QH125_3 QH125_4 QH125_7 QH125_8 QH125_10 QH125_11 QH125_12 QH125_14 QH125_16 QH125_18
  QH126_1 QH126_2 QH126_3 QH126_4 QH126_5 QH126_7 QH126_9 QH126_95 QH126_96
  QH131A
  QH132A QH132B QH132C QH132D QH132E QH132F QH132G QH132H QH132I QH132J QH132K QH132L QH132M QH132N 
  QH133A QH133C QH133D QH133E QH133F QH133G
  QH132AA
  QH133AA 
  MOBPHONE CHECKACC CREDITCARD EWALLET
  QH152_11 QH152_31 QH152_32 QH152_33 QH152_34 QH152_96
  QH153_21 QH153_22 QH153_31 QH153_32 QH153_96
  QH154_21 QH154_22 QH154_31 QH154_32 QH154_33 QH154_34 QH154_35 QH154_96
  HOUSE LAND memsleep
  /print univariate initial extraction fscore
  /criteria factors(1) iterate(25)
  /extraction pc
  /rotation norotate
  /save reg(ALL urb)
  /method=correlation.


*** Rural Factor Analysis.
filter off.
filter by rural.
execute.

* Ladys Ortiz.
* removed QH109_21_sh QH152_96.

factor /variables
  QH100A_1 QH100A_2 QH100A_3 QH100A_4 QH100A_6
  QH101_11 QH101_12 QH101_21 QH101_31 QH101_41 QH101_51 
  QH109_11 QH109_12 QH109_21 QH109_22 QH109_23 QH109_61 
  QH109_11_sh QH109_12_sh QH109_22_sh QH109_23_sh 
  QH117_1 QH117_2 QH117_3 QH117_4 QH117_95 
  QH122A 
  QH123_1 QH123_2 QH123_3 QH123_4 QH123_5 QH123_6 QH123_7 QH123_95 
  QH125_1 QH125_3 QH125_4 QH125_7 QH125_8 QH125_10 QH125_11 QH125_12 QH125_14 QH125_16 QH125_18
  QH126_1 QH126_2 QH126_3 QH126_4 QH126_5 QH126_7 QH126_9 QH126_95 QH126_96
  QH131A
  QH132A QH132B QH132C QH132D QH132E QH132F QH132G QH132H QH132I QH132J QH132K QH132L QH132M QH132N 
  QH133A QH133C QH133D QH133E QH133F QH133G
  QH132AA
  QH133AA 
  MOBPHONE CHECKACC CREDITCARD EWALLET
  QH152_11 QH152_31 QH152_32 QH152_33 QH152_34 
  QH153_21 QH153_22 QH153_31 QH153_32 QH153_96
  QH154_21 QH154_22 QH154_31 QH154_32 QH154_33 QH154_34 QH154_35 QH154_96
  HOUSE LAND memsleep
  /missing meansub
  /analysis
  QH100A_1 QH100A_2 QH100A_3 QH100A_4 QH100A_6
  QH101_11 QH101_12 QH101_21 QH101_31 QH101_41 QH101_51 
  QH109_11 QH109_12 QH109_21 QH109_22 QH109_23 QH109_61 
  QH109_11_sh QH109_12_sh QH109_22_sh QH109_23_sh 
  QH117_1 QH117_2 QH117_3 QH117_4 QH117_95 
  QH122A 
  QH123_1 QH123_2 QH123_3 QH123_4 QH123_5 QH123_6 QH123_7 QH123_95 
  QH125_1 QH125_3 QH125_4 QH125_7 QH125_8 QH125_10 QH125_11 QH125_12 QH125_14 QH125_16 QH125_18
  QH126_1 QH126_2 QH126_3 QH126_4 QH126_5 QH126_7 QH126_9 QH126_95 QH126_96
  QH131A
  QH132A QH132B QH132C QH132D QH132E QH132F QH132G QH132H QH132I QH132J QH132K QH132L QH132M QH132N 
  QH133A QH133C QH133D QH133E QH133F QH133G
  QH132AA
  QH133AA 
  MOBPHONE CHECKACC CREDITCARD EWALLET
  QH152_11 QH152_31 QH152_32 QH152_33 QH152_34
  QH153_21 QH153_22 QH153_31 QH153_32 QH153_96
  QH154_21 QH154_22 QH154_31 QH154_32 QH154_33 QH154_34 QH154_35 QH154_96
  HOUSE LAND memsleep
  /print univariate initial extraction fscore
  /criteria factors(1) iterate(25)
  /extraction pc
  /rotation norotate
  /save reg(ALL rur)
  /method=correlation.


* Label the created score variables.
rename variables (com1 urb1 rur1=comscore urbscore rurscore).
variable labels comscore 'Common wealth score' /urbscore 'Urban wealth score' /rurscore 'Rural wealth score'.

* Add a variable used for linking later.
filter off.
string rowtype_ (A8).
compute rowtype_ = 'EST'.


* Calculate regressions.
** Area=urban.
filter off.
filter by urban.
execute.

* Declare a dataset to be written to in the regression.
dataset declare urbancorv.
* Run regression of comscore with urbscore.
regression
  /missing listwise
  /statistics coeff outs R anova
  /criteria=pin(.05) pout(.10)
  /noorigin
  /dependent comscore
  /method=enter urbscore
  /outfile=corv(urbancorv).

* Activate file of output from regression.
dataset activate urbancorv.
* Drop all rows of output except the coefficients.
select if (rowtype_ = 'EST').
execute.
* Delete unnecessary variables before merging.
delete variables DEPVAR_ VARNAME_.
* Rename variables containing the constant and the coefficient.
rename variables CONST_=urbconst urbscore=urbcoeff.

* Re-activate the main household data.
dataset activate assets.
* merge the coefficients.
match files /file = *
  /table = urbancorv
  /by ROWTYPE_.
execute.


* Calculate regressions.
** Area=rural.
filter off.
filter by rural.
execute.

* Declare a dataset to be written to in the regression.
dataset declare ruralcorv.
* Run regression of comscore with rurscore.
regression
  /missing listwise
  /statistics coeff outs R anova
  /criteria=pin(.05) pout(.10)
  /noorigin
  /dependent comscore
  /method=enter rurscore
  /outfile=corv(ruralcorv).

* Activate file of output from regression.
dataset activate ruralcorv.
* Drop all rows of output except the coefficients.
select if (rowtype_ = 'EST').
execute.
* Delete unnecessary variables before merging.
delete variables DEPVAR_ VARNAME_.
* Rename variables containing the constant and the coefficient.
rename variables CONST_=rurconst rurscore=rurcoeff.

* Re-activate the main household data.
dataset activate assets.
* merge the coefficients.
match files /file = *
  /table = ruralcorv
  /by ROWTYPE_.
execute.


dataset close urbancorv.
dataset close ruralcorv.
dataset activate assets.
filter off.


*** Calculate combined wealth score from Urban and Rural Scores.
* Use coefficients from urban and rural regressions above.
compute combscor=0.
print formats combscor (f11.5).
write formats combscor (f11.5).
** Urban.
if (urban = 1) combscor=urbconst+urbcoeff*urbscore.
** Rural.
if (rural = 1) combscor=rurconst+rurcoeff*rurscore.
variable labels combscor 'Combined national wealth score'.
execute.

* Ladys Ortiz.
compute wght =  QHWEIGHT.
if (QHSYRIAN = 1 or QHSYRIAN = 2)  wght = QHSYWGHT.

*compute hhwt = QHWEIGHT/1000000.
compute hhwt = wght / 1000000.
variable labels hhwt 'HH weights'.
formats hhwt (f12.6).
weight by hhwt.

filter off.
execute.

frequencies variables=combscor comscore
  /format=notable
  /ntiles=5
  /statistics=minimum maximum stddev mean
  /histogram normal
  /order=analysis.

filter off.
filter by urban.
execute.

frequencies variables=combscor urbscore
  /format=notable
  /ntiles=5
  /statistics=minimum maximum stddev mean
  /histogram normal
  /order=analysis.

filter off.
filter by rural.
execute.

frequencies variables=combscor rurscore
  /format=notable
  /ntiles=5
  /statistics=minimum maximum stddev mean
  /histogram normal
  /order=analysis.

filter off.
execute.

* Ladys Ortiz.
*Calculate quintiles and scores for data file.
*compute hhmemwt=QHWEIGHT*hhusual/1000000.
compute hhmemwt= wght * hhusual / 1000000.
variable labels hhmemwt 'HH members weighting for index'.
formats hhmemwt (f12.6).
weight by hhmemwt.

filter off.
filter by urban.
execute.

rank variables=urbscore (A) /rank /ntiles (5) /print=yes /ties=mean.

filter off.
filter by rural.
execute.

rank variables=rurscore (A) /rank /ntiles (5) /print=yes /ties=mean.

filter off.
execute.

rank variables=combscor (A) /rank /ntiles (5) /print=yes /ties=mean.

variable labels ncombsco 'Combined wealth index' /nurbscor 'Urban wealth index' /nrurscor 'Rural wealth index'.
value labels ncombsco nurbscor nrurscor 1 'Lowest' 2 'Second' 3 'Middle' 4 'Fourth' 5 'Highest'.

*** Check on quintiles.
frequencies variables=ncombsco nurbscor nrurscor.


weight by hhwt.

ctables /table (
  QH100A_1+QH100A_2+QH100A_3+QH100A_4+QH100A_6+
  QH101_11+QH101_12+QH101_21+QH101_31+QH101_41+QH101_51+
  QH109_11+QH109_12+QH109_21+QH109_22+QH109_23+QH109_61+
  QH109_11_sh+QH109_12_sh+QH109_21_sh+QH109_22_sh+QH109_23_sh+
  QH117_1+QH117_2+QH117_3+QH117_4+QH117_95+
  QH122A+
  QH123_1+QH123_2+QH123_3+QH123_4+QH123_5+QH123_6+QH123_7+QH123_95+
  QH125_1+QH125_3+QH125_4+QH125_7+QH125_8+QH125_10+QH125_11+QH125_12+QH125_14+QH125_16+QH125_18+
  QH126_1+QH126_2+QH126_3+QH126_4+QH126_5+QH126_7+QH126_9+QH126_95+QH126_96+
  QH131A+
  QH132A+QH132B+QH132C+QH132D+QH132E+QH132F+QH132G+QH132H+QH132I+QH132J+QH132K+QH132L+QH132M+QH132N+
  QH132AA+
  QH133A+QH133C+QH133D+QH133E+QH133F+QH133G+
  QH133AA+
  MOBPHONE+CHECKACC+CREDITCARD+EWALLET+
  QH152_11+QH152_31+QH152_32+QH152_33+QH152_34+QH152_96+
  QH153_21+QH153_22+QH153_31+QH153_32+QH153_96+
  QH154_21+QH154_22+QH154_31+QH154_32+QH154_33+QH154_34+QH154_35+QH154_96+
  HOUSE+LAND+memsleep
  )[S][Mean F8.3] by (ncombsco+nurbscor+nrurscor)[C]
  /slabels visible=no.

* Alternative instead of ctables.
*means tables=
  QH100A_1 QH100A_2 QH100A_3 QH100A_4 QH100A_6
  QH101_11 QH101_12 QH101_21 QH101_31 QH101_41 QH101_51 QH101_96
  QH109_11 QH109_12 QH109_13 QH109_21 QH109_22 QH109_23 QH109_61 QH109_96
  QH109_11_sh QH109_12_sh QH109_13_sh QH109_21_sh QH109_22_sh QH109_23_sh QH109_96_sh
  QH117_1 QH117_2 QH117_3 QH117_4 QH117_95 QH117_96
  QH122A QH123_1 QH123_2 QH123_3 QH123_4 QH123_5 QH123_6 QH123_7 QH123_95 QH123_96
  QH125_1 QH125_3 QH125_4 QH125_6 QH125_7 QH125_8 QH125_10 QH125_11 QH125_12 QH125_13 QH125_14 QH125_16 QH125_18 QH125_96
  QH126_1 QH126_2 QH126_3 QH126_4 QH126_5 QH126_6 QH126_7 QH126_8 QH126_9 QH126_10 QH126_11 QH126_12 QH126_13 QH126_14 QH126_95 QH126_96
  QH131A QH132A QH132B QH132C QH132D QH132E QH132F QH132G QH132H QH132I QH132J QH132K QH132L QH132M QH132N QH132AA QH133A QH133C QH133D QH133E QH133F QH133G
  QH133AA MOBPHONE CHECKACC CREDITCARD EWALLET
  QH152_11 QH152_31 QH152_32 QH152_33 QH152_34 QH152_96
  QH153_21 QH153_22 QH153_31 QH153_32 QH153_96
  QH154_21 QH154_22 QH154_31 QH154_32 QH154_33 QH154_34 QH154_35 QH154_96
  HOUSE LAND memsleep
    by ncombsco nurbscor nrurscor
  /cells mean count stddev.


graph /histogram(normal)=combscor /title= 'Distribution of Households by Wealth Scores'.
frequencies variables=combscor
  /format=notable
  /ntiles=5
  /statistics=stddev minimum maximum semean mean median mode skewness seskew kurtosis sekurt
  /order=analysis.


weight off.
filter off.
write formats combscor urbscore rurscore (f11.5).

* Save final dataset of assets.
save outfile='assets.sav'.

*** Write out scores file.
write outfile='scores.dat' table /QHCLUST QHNUMBER combscor ncombsco urbscore nurbscor rurscore nrurscor.
execute.

* Produce table 2.6 weighted by de jure population to test.
weight by hhmemwt.
compute total = 0.
variable labels total 'Total'.
value labels total 0 ' '.
ctables /vlabels variables=QHTYPE QHREGION total Ncombsco display=both
  /table QHTYPE [C] + QHREGION [C] + total [C] BY Ncombsco [C][rowpct.count '%' F8.1, totals [rowpct.count '%' F8.1, count '' F8.0]]
  /categories variables=Ncombsco total=yes
  /titles title='Table 2.6 Wealth Quintiles' ''
    'Percent distribution of the de jure population by wealth quintiles, according to residence and region, Jordan DHS 2023'
  corner='Residence/region' /slabels visible=no.